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Abstract 



We parametrize phenomenologically possible new physics effects and calculate their 
inffuence on the invariant dilepton mass spectrum and the Dalitz plot for the decays 
B —>■ Ki~^i~ and B —>■ K*i~^i~. Especially the decay into K* yields a wealth of new 
information on the form of the new interactions since the Dalitz plot is sensitive to 
subtle interference effects. We also show how transversely polarized i^*-mesons give 
additional information. 
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1 Introduction 



Rare decays of _B-mesons, such as the recently observed processes B K*"-/ and 
B Xs'~f 0] may become an important tool for studying new forces beyond the 
standard model. Indeed, many authors have investigated the effects of multi-Higgs 
models, super symmetric theories, left-right symmetric models, etc. on this process 
. The interest in these decays stems from the fact that they occur in the standard 
model only through loops and are therefore particularly sensitive to "new physics" . 

Whereas most of the past work has been concernend with the decay b ^ s'y, the 
leptonic transition b — > si~^i~ is more sensitive to the actual form of the new inter- 
actions since it not only allows to measure a total rate, but also various kinematical 
distributions. More specifically, while new forces will only change the normalization 
of 6 — > S7 (or B K*'~f), the interplay of the various operators will also change the 
spectra in the decays with two leptons. Thus, different new physics models which 
yield the same 6 — >■ 57 amplitude may be distinguished by their characteristic invari- 
ant mass spectrum of the lepton pair in 6 — ^ s£~^£~ . Of course, there is a spectrum 
for the photon in the inclusive process B X^^j, but it is essentially due to the 
wave function of the decaying meson and therefore not of interest here. 

In this note we will parametrize possible new physics in a phenomenological way 
and evaluate their effects on the exclusive decays B Ki~^£~ and B — > K*i~^i~. 
Similar results also follow for the inclusive lepton pair spectrum, but since coherence 
effects tend to be more pronounced in exclusive decays we have concentrated here 
on these. Of course, this introduces a certain dependence on the hadronic model 
used. However, since most models do not disagree violently in the regions of in- 
terest and the crucial signs between the form factors are quite safe, we feel that 
a reliable representation of the effects is indeed possible. We have made use of a 
recent quark model 0, ||, P] which gives rather satisfactory results. Needless to say 
that it reproduces the heavy quark limit. 



2 Effective Hamiltonian 

The effects of both the new and the heavy standard model particles can be summa- 
rized by an effective Hamiltonian; for b — ^ si~^i~ (and the corresponding exclusive 
decays) it has the form Q 

^ = ^E^^(/^)o.(/i) , (1) 

where Gp denotes the Fermi constant and Xt = VtbV*^ are products of CKM matrix 
elements. The Ci{fi) are coefficients calculated from processes with heavy particles 
exchanged and Oi are the following operators 

01 = (Sa'J^Lba) {cp'^^Lcp) , 

02 = {Sal^Lbp) (CfS'Jf.LCa) , 
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g 



(2) 



The basis of ref. has been supplemented by the primed operators; the operators 
not written give small contributions. In eq. (|1]), the renormalization scale n is usually 
chosen to be /x ~ rub in order to avoid large logarithms in the matrix elements of 
the operators Oj. Accordingly, the Ci are shifted from their perturbative value at 
/i ~ mw by the renormalization group equations. As the analytic expressions for the 
coefficients Cj(/i) in the standard model are given in many places in the literature 
0, we only note here their numerical values which correspond to = 174 GeV 
||TD|: At the scale /i = mw one gets 



, C2{mw) = -1 , Crimw) = 0.196 



Ci{mw) 

171, ^ 

, [mw) , Cs{mw) = -2.02 

rrib 

Cg{mw) = 4.56 , Cs'{mw) = C9>{mw) = , 
while at /i = the non- vanishing coefficients are 

Ci(mb) = 0.226 , Ciirrib) = -1.096 , Crirrib) = 0.310 



(3) 



Crinib) 



nib 



Cj{mb) , CsK) = -3.84 , C^{mb) = 4.56 



(4) 



It is well-known that the operators Oi and O2 generate both short and long distance 
contributions to 6 — » si'^i" and to the corresponding exclusive decays. The long 
distance contributions are mainly due to the J/ip and ip' resonances jn], |T2|, 
These short and long distance contributions can be easily built-in [|, by replacing 
the coefficient by Cg^^ , i.e. 



Cl^'irub) = C8(m,) + [3Ci(m,)+C2K)] X 



h{rhc, s) - 



a 



E 



em Vi=J/il),^' 



(5) 



where = rric/mb and s = /ml with q^ beeing the invariant mass squared of the 
lepton pair. The short distance contributions are contained in the function h{rhc, s) 
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which reflects the one- loop matrix element of the four-quark operators Oi and O2. 
For s > Az'^ 



h{z, s) 



4, . 8 16^2 2 / 4^2 / 4^2N 

-logz^ h 2 H 

9^ 27 9 s 9V si s, 



X 



l0£ 



1 + ^1 - Az-^/s 



1 - Jl-Az^/s 



+ in 



(6) 



and for s < Az^ 



h(z, s) 



log 



8 16 z^ 
27 ~ YT 



■+ 



4 4^2 / 4^2 \ 

2 + - atan 



/4^2/s - 1^ 



(7) 



The value of k in eq. (^) must be chosen such that the combination 

fi:[3Ci(mb) + C2K)] ^-1 

in order to correctly reproduce the branching ratio [|l4| 

BR{B -> J/^X trX) = BR{B -> J/^X) BR{J/ip tt) . 

In this language, the physics beyond the standard model is characterized by the 
values of the Wilson coefficients Ci at the perturbative scale my/- For instance, 



supersymmetry can lead to a C^ which is opposite to its standard model value [15 



whereas left-right symmetric models could give a range for C7/ between minus and 
plus the standard model value of C7 0, |I6|, |T^. Models with leptoquarks may cause 
the other coefficients to lie within similar intervals. In order to limit ourselves to a 
reasonable set of values in the examples, we chose to vary the coefficients between 
minus and plus their standard model values. In addition, we require that the rates 
for B — ^ and B K*^ do not change. For both reactions, this implies that 
we must leave the value of (|C7(m;,)|2 -|- |C7'(mfe)|2) fixed within certain limits. 

The running of the coefficients C7/, Cg' and Cg/ is simpler than the one of the 
unprimed coefficients since there is no operator mixing with O2. In principle, new 
operators Oi' and O2' with right-handed four fermion interactions should be in- 
cluded; we assume them to be small, as indicated by the high accuracy of the V — A 
form in many weak decays, and we neglect them. For the scaling we can simply 
write in leading-log accuracy 



Cg, (mfe) 
Cr {nib) 

where i] = as{mb) / as{mw) — 1-71 for nib = 5 GeV. 



Cs'imw) 



(8) 
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3 The decay B K£+£ 

The relevant hadronic matrix elements for the decay B Ki~^i~ associated with 
the operators Oi are parametrized as in eqs. (2.11) and (2.12) of ref. [Q, i.e. 



{pK\siaf,uq''{lTl5)b\PB) 



1 



(9) 
(10) 



mB + rriK 

where P = Pb + Pk and q = Pb — Pk and 0"^,^ is defined as a^u = {i/'^)[y^l, lA- The 
form factors F+, F_ and Ft are written generically as 



F{q' 



F(0) 



(11) 



The parameters Ai and A2 characterizing the shape of each of the form factors are 
given together with -F(O) in table 1 which has been taken from ref. f^. The full 
information on the decay depends on two kinematical variables which we choose to 
be s = q^ /ml and x = E^/ttib, where is the invariant mass of the lepton pair and 
Ei is the energy of the negatively charged lepton measured in the rest frame of the 
B meson. Using the notation of ref. [0 we get the double differential decay width 



dxds 



Glm%\Xt 
87r3 



2 

~ 167r2 



(3 [4:X{Xm - x) - s(l - 2x)] 



where 



^B + rriK 

The variables x and s vary in the range 



+ |(C9 + C90i^+ 



12) 



(13) 



4m^ „ {rriB 



rriK] 



s + 2x„ 



70 . S + 2X^ + a/0 

< X < ■ 



B 



nil 



with 



Xv] 



2m| 



is + 2x^Y - As. 



(14) 



(15) 



Integrating eq. (0) over the variable x in the range specified in eq. (0) yields the 
invariant mass distribution of the lepton pair 



dr{B^Ki+i~) Glm%\Xt\^ 



a: 



ds 



(16) 



We notice the relative sign in eq. (0) between the terms involving the Cj and 
Cg coefficients which leads to a potentially interesting interference effect. Since -F+ 
and Ft are calculated in a specific model, one might worry that the sign is strongly 
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model dependent. However, in the approximation where the quarks are on-shell, the 
Gordon decomposition confirms the sign. 

In Fig. 1 we show the dilepton invariant mass distribution for the standard model 
and for a model where Cg has the opposite sign from the standard model value. The 
difference between the two curves is especially noticeable near the resonances, a 
fact which is easily explained by eq. (^) (see also eq. (|T3p). On the other hand, 
changing the sign of Cy does not alter the standard model picture very much, since 
the contributions from the C7 terms in eq. ([TB|) are small. From the form of (3 we 
see that other variations of the coefficients will not give dramatic changes of the 
spectra and Dalitz plot. 

In the numerical evaluations we take nib = = 5.28 GeV in all formulas except 
in the Wilson coefficients, where we used = 5 GeV. Furthermore, in all the plots 
we divided the spectra by the total i?-meson decay width, estimated to be [|1^ 



tot 



fmlGl\V,b\ 
1927r3 



(17) 



with / ^ 3. 



4 The decay B K*i+t 



Next, we consider the more interesting decay B K*i'^i~. Again, the relevant 
hadronic matrix elements are parametrized as in eqs. (2.18) and (2.19) of ref. [^, 
i.e. 



{PK*\s'J^,{lTl5)^PB) 



1 



rriB + rriK 

±Ao(g2)(m| - m|0< ± A4q')ie*P)P, ± A4q')ie*P)q, 



(18) 



{pKAs^(r^,q%l ± -f^Ws) = -ig{q^)e^uaf3e*'' P'^q" ± 



ao(g^ 



B 



m 



Kj 



±a^{q') (e*P) 



(19) 



where P and q are now defined as P = + px* and q = Pb — Pk*- The e- 
tensor is taken in the Bjorken-Drell convention (i.e. eoi23 = 1) contrary to ref. 0, 
which explains the additional minus sign of the terms proportional to the e-tensor. 
The various form factors in eqs. (0) and (p!9| ) are written as in eq. (0) and the 
parameters are given in table 1. With the notation of ref. |19| we get 



K* 
dxds 



Glm%\Xt\'^ qL 
-27s(xm - 



— a + 2l3[Ax{xr. 



m 



B 



2x + s/2)] 



x) - s(l - 2x)] 

(20) 



6 



where 



a 
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\G'a\' + \n\' + ^rnlpl {\Gy\' + \Fy 



\G'a\' + \F1\ 



2 

ml 



- mis {\Gy\^ + \Fy\^) + {\G\^ + \Ft\') + 



K* 



-2Re(G° Fy + F\G 



V 



(21) 



and 



Z J, -O 

Ps=<P— ■ 



The functions Gy, Fy, G% F^, Gj^ and F^ are combinations of form factors and 
Wilson coefficents; they read 



Gy 
Fy 



Fl 



G1 



2(m. + m.) 



2{mB + m/^.) 



(Gg — Gg 



I 



2{mB + mK*) q 
(Gg - G9OA+ 



2(mB + rriK* 



(22) 



The range of the variables x and s and the definition of obtained from eqs. 

(IT^ ) and (15) by the obvious replacement mx —>■ mx*- 

Integration over the variable x leads to the invariant mass distribution of the 
lepton pair 



dr{B -> K*l 
ds 



167r2 



a -\ — 



m. 



(23) 



Taking into account only transversely polarized K* mesons in the final state, the 
analogous expression for the double differential decay width is 



d^T^{B K*i+i-) 



dxds 



47]-3 



167r2 



m. 



■ a' + 2/3' [Ax{x 



X 



-s(l-2x)] 



-27^s(x^ - 2a; + s/2) - 



2ml' 



'4:x{xm — x) — s(l — 2a;)' 

, {Xm + s/2f - S ^ 



, (24) 



where 



a 



a 



0" 



-m 



B 



s{\Gv\^ + \Fv 



7 



7 



?0\2 



(25) 



The corresponding invariant mass distribution of the lepton pair reads 



dV^iB K*i+t 
ds 



Glm%\Xt\^ 



12m| 



at 



167r2 



■a + — — 



mi 



(26) 



In Figs. 2 to 4 we show the invariant mass distribution and the Dahtz plots for 
three sets of values of the couplings, all of which give the same rate for the decay 
B —>■ Xsj. In Fig. 2 we illustrate the standard model prediction (Fig. 2a) and 
the case where the sign of Cg is reversed (Fig. 2b); this again leads to the marked 
deviation around the resonances. In addition, we see from eqs. (^Tj) and (^3]) that 
reversing the sign of Cg does not change the invariant mass distribution; however, 
this change influences the Dalitz plot significantly as illustrated in Fig 3: While the 
standard model distribution in Fig. 3a is populated more for large x-values. Fig. 3c 
shows more points for small values of x. 

Finally, from Fig. 4 we see that the transverely polarized K* mesons accentuate 
this asymmetry in the Dalitz plot and can therefore be used to fix the sign of the 
coefficient Cg. 



5 Conclusions 

In this paper we have investigated the effects of new physics on the rare decays 
B — ^ Ki~^i~ and B K*i~^i~. Whereas the new physics influence only the rate 
in the decay B Xg'y, they change signiflcantly the various spectra in the leptonic 
decays. Thus, while certain models give the same rate for B Xg'y, they can be 
distinguished by their dilepton spectra. In addition, the Dalitz plot carries further 
information (compare Figs. 3 and 4). Furthermore, looking only at the transversely 
polarized K* yields improved information; as an example we have investigated the 
coefficient Cg. 

Of course, the applicability of our results will depend on the number of available 
decays. It is expected that at the hadronic B facilities at LHC the required number 
(about 100 for the dilepton spectrum and 1000 for the Dalitz plot) can be analyzed. 
Then, these decays will be a truly powerful tool towards understanding new physics. 



Noted added: A. Ali, G. Giudice and T. Mannel have investigated in a similar 
spirit the inclusive dilepton decays (private communication and ||20|| ) . 
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Figure captions 

Fig. 1: Invariant dilepton mass distribution for the decay B K£^£^ . In Fig. la 
the prediction within the standard model is given; Fig. lb shows the result 
for a model where Csimw) = 2.02 (opposite value of the standard model; see 
eq. (§)) and the other coefficients unchanged. 

Fig. 2: Same as Fig. 1, but for the decay B K*£~^i~. The solid (dashed) line 
takes into account all (only the transverse) polarizations of the i^'*-meson. 

Fig. 3: Dalitz plot (compare eq. (PP])) for the decay B K*i^i^ taking into 
account all polarizations of the i^*-meson. We have used 8000 points in the 
scatter plots. The distribution in Fig. 3a is for the standard model values of 
the coefficients while Fig. 3b has the sign of Cs{mw) reversed. In Fig. 3c we 
used the standard model values of the the coefficients, except that the sign of 
Cg^mw) has been reversed. 

Fig. 4: Dalitz plot for transversely polarized i^*-mesons. Fig. 4a is for the standard 
model values of the coefficients and Fig. 4b has the sign of Cg{mw) reversed. 

Tables 





F+ 




V 




^0 


9 


a+ 


ao 


m 


0.30 


-0.30 


-0.35 


0.24 


-0.37 


0.31 


-0.31 


0.31 


Ai (GeV) 


4.06 


4.08 


4.05 


4.39 


5.92 


4.05 


4.37 


6.40 


A2 (GeV) 


5.51 


5.50 


5.45 


5.81 


8.54 


5.44 


5.84 


8.84 



Table 1: Values of the various F{0), Ai and A2 
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